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We present a theoretical model that describes electrical spin-detection at a ferromag- 
net/semiconductor interface. We show that the sensitivity of the spin detector has strong bias 
dependence which, in the general case, is dramatically different from that of the tunneling current 
spin polarization. We show that this bieis dependence originates from two distinct physical mecha- 
nisms: 1) the bias dependence of tunneling current spin polarization, which is of microscopic origin 
and depends on the specific properties of the interface, and 2) the macroscopic electron spin trans- 
port properties in the semiconductor. Numerical results show that the magnitude of the voltage 
signal can be tuned over a wide range from the second effect which suggests a universal method for 
enhancing electrical spin-detection sensitivity in ferromagnet/semiconductor tunnel contacts. Using 
first-principles calculations we examine the particular case of a Fe/GaAs Schottky tunnel barrier 
and find very good agreement with experiment. We also predict the bias dependence of the voltage 
signal for a Fe/MgO/GaAs tunnel structure spin detector. 

PACS numbers: 



I. INTRODUCTION 

Semiconductor spintronics aims to harness the elec- 
tron's spin degree of freedom in data storage and process- 
ing, typicaUy by utihzing heterostructures composed of 
a combination of magnetic and non-magnetic materials^. 
A fundamental problem in semiconductor spintronics was 
to find ways to electrically generate non-equilibrium elec- 
tron spin distributions in conventional semiconductors. 
Efficient electrical spin injection from ferromagnetic con- 
tacts into semiconductors using spin dependent tunnel- 
ing ^, ^ was shown to overcome the 'conductivity mis- 
match problem' associated with highly conductive metal- 
lic contacts^. Spin polarization of the tunneling current 
originates from the spin dependence of the electron wave- 
functions and the densities of states of the ferromagnetic 
contact. The spin dependent tunneling approach was 
realized experimentally using Fe interfaces with GaAs 
^j6j7j8j9^ siliconiSiiii and graphenSiiS.. Jiang et. al using 
CoFe/MgO interfaces showed enhanced electrical spin in- 
jection efficiency into GaAsH. Contacts made of CoFe 
were also used to inject directly into GaAsii and showed 
an electron polarization that had dramatically different 
bias dependence from that of Fe/GaAs contacts^. In ad- 
dition to electrical spin injection, efficient electrical spin 
detection is required to achieve functional semiconductor 
spintronic devices. 

Crooker et al^ recently reported experiments of elec- 
trical spin detection using Fe/GaAs Schottky tunnel bar- 
riers as electrical spin detectors. They demonstrated that 
both the magnitude and sign of the spin detection sen- 
sitivity are tunable with voltage bias applied across the 
Fe/GaAs interface; in some cases they were able to im- 
prove the spin detection sensitivity by an order of mag- 
nitude. The bias dependence of the sensitivity of the de- 
tector was shown to be dramatically different from that 
of the injected current spin polarization. A theoretical 
model was used to correlate the spin-detection sensitiv- 



ity of the Fe/GaAs electrodes with their bias-dependent 
spin injection properties. The model described success- 
fully many of the experimentally observed trends^^. 

Here we give a detailed description of this theoret- 
ical model of electrical spin detection at a ferromag- 
net/semiconductor interface. We consider a case when 
spin polarization is generated in the semiconductor by 
an external source (e.g., optical or electrical) and subse- 
quently detected at a ferromagnetic contact in which a 
tunnel barrier exists at the ferromagnet/semiconductor 
interface. We incorporate first principles calculations to 
examine two specific cases of tunnel barrier, a Fe/GaAs 
Schottky tunnel barrier and a Fe/MgO/GaAs tunnel 
structure. In this way we demonstrate that the theory is 
general and can be applied to a variety of electrical spin 
detectors. We show that the sensitivity of the electrical 
spin detectors have strong bias dependence of both mi- 
croscopic and macroscopic origin. While the bias depen- 
dence of microscopic origin is specific to each ferromag- 
net/semiconductor tunneling structure, the macroscopic 
bias dependence is general and depends on the electrical 
transport properties of the semiconductor. 

This article is organized as follows. In Section II we 
give a detailed description of the theory. At the end of 
the section we provide some approximate analytical for- 
mulas and look at the asymptotic limit of large currents. 
This helps in understanding the general trends predicted 
by the theory and facilitates a discussion of the numerical 
results. In section III, first we present numerical results 
for a generic spin detector with a tunneling current polar- 
ization that has no bias dependence, then we incorporate 
first principles calculations to examine the specific cases 
of Fe/GaAs Schottky barrier and Fe/MgO/GaAs tunnel 
structure spin detectors. We conclude the paper in Sec- 
tion IV. with a summary of our results and conclusions. 
Some calculational details are included in an appendix. 
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II. THEORETICAL MODEL 

A voltage signal, in response to a change in the spin po- 
larization of the current, at a ferromagnetic tunnel junc- 
tion results because the tunneling resistance of the junc- 
tion depends on electron spin. For example, if the tun- 
neling resistance of junction is smaller for spin-up (ma- 
jority spin in the ferromagnetic contact) electrons than 
for spin-down electrons (minority spin in the ferromag- 
netic contact), the voltage drop across the junction will 
decrease if the spin-up electron component of the cur- 
rent increases while the spin down electron component of 
the current decreases keeping the total current constant. 
In recent experiments reported by Crooker et al^^, the 
spin polarization of the current at the tunnel junction is 
changed either by the absorption of circularly polarized 
light in the semiconductor in the vicinity of the ferro- 
magnetic tunnel junction or by electrical injection from 
a remote ferromagnetic contact. The spin polarized elec- 
tron density generated by the absorption of the circularly 
polarized light or remote electrical injection drift/diffuses 
to the detection ferromagnetic tunnel junction and thus 
modifies the spin polarization of the current across this 
junction. A bias dependence of the voltage signal, as is 
observed experimentally, occurs because of a combina- 
tion of two effects: 1) transport of the remotely gener- 
ated spin polarized electron density depend on bias so 
that the change in spin polarization of the current cross- 
ing the tunnel junction depends on voltage bias; and 2) 
the spin dependence of the tunneling resistance depends 
on voltage bias. 

Fig. [1] shows a schematic diagram of the electro- 
chemical potentials for spin-up and spin-down electrons 
in the vicinity of a ferromagnetic tunnel junction in which 
the junction resistance is smaller for spin-up electrons 
than for spin-down electrons. We show for the cases un- 
der reverse bias (electron injection into the GaAs) and 
under forward bias (electron extraction from the GaAs), 
without remote spin generation (solid lines) and with re- 
mote generation of spin- up electrons (dashed lines). The 
Fe contact is highly conductive so that the electric field in 
the Fe is very small and the spin-up and spin-down elec- 
tron distributions in the Fe are very nearly in equilibrium 
with each other. Thus, in the Fe, the electro-chemical 
potentials for the two spin types are nearly constant in 
position and equal to each other independent of bias. 
The shaded stripe in Fig[T] represents the depletion re- 
gion in the GaAs that forms the tunnel barrier. There is 
a drop in the electro-chemical potential across this tunnel 
barrier region for both spin-up and spin-down electrons. 
This drop in electro-chemical potential is proportional to 
the product of the current and the tunnel resistance for 
each spin type. It is larger for spin-down electrons than 
for spin-up electrons because the tunneling resistance is 
larger for spin-down electrons than for spin-up electrons. 
In reverse bias, the electro-chemical potentials are higher 
in the Fe contact than in the GaAs and because the drop 
across the junction is larger for spin-down electrons than 
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FIG. 1: Schematic of near a Fe/GaAs Schottky de- 
tector, at forward and reverse bias. Dotted lines show how 
fil^ l are modified by an additional, remote source of spin po- 
larization. In reverse bias remote spin generation increases 
— Mi) near the contact but not as much as in the region 
of spin generation. 



for spin-up electrons there is a surplus of spin-up elec- 
trons compared to spin down electrons near the junction 
interface. In forward bias the electro-chemical potentials 
are higher in the GaAs than is the Fe contact and there 
is a surplus of spin-down electrons compared to spin-up 
electrons near the junction interface. Because of spin re- 
laxation, the electro-chemical potentials for spin-up and 
spin-down electrons come together far from the tunnel 
junction when there is no spin generation. When there is 
generation of spin-up electrons, the two electro-chemical 
potentials far from the tunnel junction remain separated 
and are determined by the spin generation profile. At 
any point in space the spin polarization of the density 
is determined by statistics from the electro-chemical po- 
tentials and the spin polarization of the current is de- 
termined by the product of the density and the spatial 
derivative of the electro-chemical potential for each spin 
type. The voltage signal at the tunnel junction produced 
by optical spin generation depends on the change in spin 
polarization of the current at the tunnel junction and the 
tunneling resistance of each spin type. 

To describe theoretically the voltage signal wc use the 
1-dimensional model of Rcf. 3 with some modifications. 
In this model the current fiow at the interface is described 
using a spin dependent interface conductance: 



(1) 



where is the current density at the interface, is the 
interface conductance, A/i,, the interfacial discontinuity 
in electrochemical potential for electrons with spin pro- 
jection ?7, and e is the magnitude of the electron charge. 
We set the interface at the position x — with the fer- 
romagnet on the left side {x < 0) and the semiconductor 
on the right (a; > 0). As in Ref. |^, we define a variable 
f3 such that = Pj, where j is the total current density 
(independent of position in the 1-dimensional model) and 
a variable a such that n-f = an, where n = n| + n| is 
the electron density in the semiconductor which is inde- 
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pendent of position and of bias. Eq. ([T]) can be written 
as: 



(A.f(0)-A.f(0)) = ^ 



(2) 



and 



(/if(0)-Mf'(0)) 



G, ' 



(3) 



where iJ,^'^{0) is the electro-chemical potential for elec- 
trons at the right (left) side of the interface and is /3 
evaluated at the interface. Adding equations ([2]) and ([3|) 



gives 



(A*f(0) + /if(0)) = (A*t(0)+A^i(0))- 

and subtracting them yields 

(A.f(0)-Mf(0)) = (A*t(0)-Mi(0)) + 
1 1 \ 1 



(4) 



The voltage drop at the interface is 



(5) 



V = [(^f (0) -h (0)) - K(0) + Mf (0))] /2e (6) 
so that 



V={j/2){(i' 



Gj G 



1 

gT 



(7) 



We consider a case in which the electrons in the semi- 
conductor are strongly degenerate and relate the electro- 
chemical potentials and the densities using a zero tem- 
perature Fermi function, then 



Aif (0) = 2"^Ef (a°) 



2/3 



and 



(0) = 2"''Ef (1 - a") 



0n2/3 



(8) 



(9) 



where Ep is the Fermi energy and a'^ is a evaluated at 
the interface. 

We describe the current flow in the semiconductor us- 
ing spin dependent drift-diffusion equations 



j-n = ^^~g^ + e^iErijj 



(10) 



where D is the diffusion coefficient, E is the electric field 
in the semiconductor (independent of position in the 1- 
dimensional model), and n is the electron mobility. The 
drift-diffusion equation at the interface gives: 



e^EAn{0) (11) 



If electrons with different spins are driven out of local 
quasi-thermal equilibrium at some region in space, so 
that rt| is not equal to n^, the difference in the two elec- 
tron densities An(x) — (n^ — rtj^)(x) relax as described 
by the spin current continuity equation 



D 



d^An 
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dAn 

dx 



f{x) 

Ts 



(12) 



where is the spin-relaxation time and f{x) is the spin 
generation function. 

We are interested in how the interface voltage drop 
varies with small changes in the amplitude G of the spin 
generation function f{x) = [C /ts)F{x) where F{x) is a 
unitless and normalized function of position. The voltage 
signal depends on the amplitude of the spin generation 
function through the spin polarization of the current at 
the interface Ref. M 



— |c=o~(l/2) [g-^-gI)J^\c=^ 



(13) 



The voltage signal depends on the derivative of 
^^^\x=o and An(0) with respect to C and to calcu- 
late these derivatives we use Eq. p^ . We search 
for a solution of this equation of the form An(x) = 
Q{x) + f{x')g{x,x')dx' subject to boundary condi- 
tions An(0) = Qo and An(oo) = where Q{x) = 
Qoe~^^^+ is the solution of the homogeneous differen- 
tial equation and 17(2;, x') is the Green's function solution 
of 
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d'^g{x,x') dg{x,x') g{x,x') 



dx"^ 



fiE 



dx 



= -S{x-x') 



(14) 

with boundary conditions ^(O, x') — and 5(00, x') — 0. 
For a general spin generation function, we have 



dAn. 
dx 

where 



An{0) 



/(^')^^^U=odx' (15) 
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, x' < x 



(16) 
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and 
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(17) 

Using j = (jE, where a = e^m is the conductivity of the 
semiconductor, As ~ y/Drl and the Fermi hquid rela- 
tionship D/fi = 2Ep/3e, we can write A± as 



1 



1 



A±(j) A, 



/3ejA, ^ 



\AaEF J 



(18) 



To be specific, we consider a striped spin generation func- 
tion of the form: 



^ ' 10 otlierwise ^ ' 

where y is the center position of striped spin generation 
and d is its width {y > d/2). Then Eq. (fTT|) can be 
written as 



of bias. Eq. (j23p is non-hnear and is solved numerically 
for a° after elimination of with the help of Eq. (|22p . 
The spin-detection sensitivity, defined as dV/C, and the 
current polarization, 2/3*^ — 1 , are then calculated numer- 
ically. 

Before we present the numerical results of the model, 
it is instructive to linearize Eq. (j23p . in order to obtain 
an analytic expression for dV/dC valid for weak injection 
conditions. For the linear case it is convenient to rewrite 
Eqs. ([21]) and in terms of a parameter 5° = Q!° - 1/2 
because a° is close to 1/2: 



(2/3° - l)j 



AEpcr 
3eA+ 



[^0 ^ gdp^^^^ ^ 250j (24) 
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where (^"^ = We eliminate /?" from the second equa- 
tion and differentiate with respect to Sd 



nO _ „-0 



eD 



[An{Q) + CF{j)] + eidEAn{Q) (20) 



where in general 



' dS^ 



2EFa 
' 3eA+ 



F{x)e''^-d{x/k^). 



(21) 



For the striped spin generation function F{i) = 
gy/A- (^grf/aA- _ g-d/2A_~j g^j^j depends on current 

through A_. (F(j) is negative because A_ is negative). 

Because + = e^iEn, - = (2/3" - \)j and 

An(0) = n(2Q;'' — 1) this can be written in terms of the 

values of parameters a and (3 at the interface. 



(2/3" - l)j = 
2EFa 



3eA 
(2a" - l)j 



[{2a°-l) + {C/n)F{j)] 



(22) 



Considering that on the metal side of the interface the 
electrochemical potentials for spin "f and | electrons are 
very nearly equal, Eqs. (O, ([5]) and (O give 



2^Ef 
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Gt G 
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(23) 



Equations and ([231) together give /3" and a" for any 
given value of the amplitude G of the source function and 



d6°_ 
dS'^ 



(27) 



where 



2t 



(gt + g\) 



(28) 



When the density polarization is small (5" ^ and 



1/2 + (5")"'/^+ (1/2 -,5") 
we can linearize A(j) 



-1/3- 



1 



2'*/3. In this case. 



(29) 
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The current polarization is given by, 



(2/3" - 1) 



A{j) 



(30) 
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We note in passing that equations (|30p and ([7]) combined 
can be used to extract Gt and G\ from the exerimental 



data set of total current, voltage and current polarization. 

Equations ((29l) . ((27|) and ((26|) give an analytic expres- 
sion for dV/dC, which in terms of the parameter 5'^ can 
be written as 



so that 



d6° _ F(j)R+ 3ej 
dS'i 4 2Ef 



(35) 



lc=o = (l/2n)^|,.=o = (l/2n) ( — - — ) 



Gt G 



(31) 

In the derivation of Eq. (|26|) we kept the diffusion 
(first RHS term) and drift (second RHS term) contribu- 
tions to the voltage change separate. The diffusion term 
has two contributions. The term —{2EF<j/3eA^)F{j) 
is always positive. It describes the change of current 
polarization at the interface, due to diffusion, resulting 
from the Green's function term in Eq. (jlSp . The term 



— (2£'Fcr/3eA+)^ is always negative because ^ is pos- 
itive. It describes the change of current polarization at 
the interface, due to diffusion, resulting from the homo- 
geneous solution term in Eq. ([T5]) The drift term, jfj^, 
describes the change in spin current polarization at the 
contact from drift of the modified spin density polariza- 
tion at the interface due to external spin generation. Its 
sign depends on the sign of j . For spin injection (j < 0) 
the absolute value of ^ cannot be larger than |F(j)| 
and hence the diffusion term is always positive. The drift 
term is always negative and therefore the two processes 
always oppose each other for spin injection. For spin 
collection (j > 0) ^ is smaller than |F(j)| for small 
currents but can become larger than |-F(j)| for large cur- 
rents. As a result the diffusion term is positive for small 
currents, but can become negative for large currents and 
the drift and diffusion augment each other for small cur- 
rents but can oppose each other for large currents in spin 
collection. 

It is interesting to examine the asymptotic behavior 
for large biasii of dV/dC. When 3e |j| As/i4:aEF) > 1, 
we can write 



A+(j) 

For the case of spin collection, where j is positive, this 
becomes as 



3ejA, ^ 3e|j|A, , 1 



o 3e|j|A 



(32) 



1 3ej 
A+(j) ^ 2aEF' 



(33) 



Then we have 



where R+ = (1/G| + 1/Gx) and 



(34) 



■' d6d 8Ef 



R+j] (36) 



Part of the diffusion term (RHS first term) cancels ex- 
actly the drift term (RHS second term) leaving a linear 
dependence on j, 



^^|,.„.-(V2„)(i--i-)F(i)i (37) 

Note that F(j) is always negative and for the case of a 
homogeneous source F(j) = —1. 

For the case of spin injection, where j is negative, Eq. 
(l32t can be written as 



1 _ 1 2crEF 
A+(j) ATSe |j|A, 



Then, 



and 



3e|j| 4 2EFa^ 
2E^^Rl^3e \j\Kl 



(38) 



(39) 
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'WA{j)ieKl IjI 



(40) 



so that 
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3e|j| 
2Ef 



(41) 



The first two RHS terms are the diffusion and the third 
the drift contributions to the change of current spin po- 
larization at the interface due to external spin genera- 
tion. The leading order terms in \j\ have opposite signs 
and cancel. Then to leading order in |j|. 



.d^ 

^ d5<^ 
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(42) 
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FIG. 2: (Color online) The calculated spin detection sensitiv- 
ity for an electrical spin detector with various levels of doping 
and tunneling current spin polarization (solid red/gray line) 
which is independent of bias. The vertical axis is logarithmic. 
The left panel is for /i — 3000cm^V~^ and the right for 
fi = 1500cm^ V~^s~^. The values of carrier concentrations 
used are depicted in the plot. The spin-relaxation time, Ts, is 
set to 10~^s. The source is set at a distance of 2 x 10~^cm 
from the interface and it has a width of 5 x 10~'*cm. The 
source function is constant within this interval. 



FIG. 3: (Color online) The calculated spin detection sen- 
sitivity for an electrical spin detector with tunneling cur- 
rent spin polarization which is independent of bias (solid 
red/gray line). The vertical axis is logarithmic For the case of 
H = 3000cm^y"^s"^ andn = 1 xlO^'^ cm"^. We vary the tun- 
neling resistance by multiples of 2. The spin-relaxation time is 
set = 10"^s. The source is set at a distance of 2 x 10 cm 
from the interface and it has a width of 5 x 10~''cm. The 
source function is constant within this interval. 



and 



dV 
dC 



|c=o ~ -(l/2n) 



1 

4 



(43) 



Unlike spin collection where the spin detection sensitivity 
grows linearly with during spin injection spin sensi- 
tivity drops as 1/ I j I . It is worth mentioning that the ex- 
act cancellation of leading order terms in Eq. (|4ip occurs 
strictly in the ID case, but there is no physical principle 
that demands exact cancellation for a more complicated 
geometry. This may lead to a reversal of the sign of the 
voltage signal with bias for spin injection. 



III. RESULTS 

We present results of the model from numerically solv- 
ing the nonlinear equations. In the following, in order to 
facilitate the comparison with experiments we adopt the 
convension that negative voltage corresponds to spin col- 
lection and positive voltage to spin injection. Often, as is 
done in experimental works we will refer to the negative 
voltage as forward bias and to the positive voltage as re- 
verse bias. In Fig.[2]we show the calculated spin detection 
sensitivity for an electrical spin detector with an interface 



which has a constant current spin polarization (indepen- 
dent of bias). Required inputs to the calculation of the 
voltage signal are: electron density, mobility, spin life- 
time, and the spin tunneling conductances G|,|. The left 
panel in Fig. [2] is for mobility /i = 3000 cm^V~^s~^ and 
the right for fj, = 1500 cm^V~^s^^ . We have used several 
values for carrier concentrations and they are shown in 
the plot. The spin-relaxation time is set to = 10~^s. 
The source is set at a distance of 2 x 10^"^ cm from the 
Fe/GaAs interface and it has a width of 5 x 10~'*cm, the 
source function is constant within this interval. As seen 
in Fig [21 in all cases, for spin injection the magnitude of 
the calculated voltage signals drops rapidly with increas- 
ing bias even though the current polarization is constant. 
By contrast, the magnitude of the calculated voltage sig- 
nals, increases rapidly with increasing bias during spin 
collection. The difference between the calculated voltage 
signal and current polarization is larger for smaller val- 
ues of the electron mobility and electron concentration. 
The magnitude of the voltage signal is smaller than the 
current polarization in reverse bias but larger in forward 
bias for a combination of two reasons: 1) the drift and 
diffusion contributions to Eq. ((26)l oppose each other 
in spin injection but add in spin collection; and 2) the 
electric field in the semiconductor tends to drift the op- 
tically or electrically generated spin polarized electrons 
away from the detector contact in reverse bias but toward 
the detector contact in forward bias. When the semicon- 
ductor is more heavily doped the electric field is smaller 
therefore these effects are less pronounced than when the 
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FIG. 4: (Color online) The calculated spin detection sensi- 
tivity for an Fe/GaAs Schottky electrical spin detector with 
various levels of doping. On both panels, the red/gray solid 
line with circles is the Fe/GaAs tunneling current spin polar- 
ization calculated from first-principles. The left panel is for 
/i = 3000cmV~^s~^ and the right for fi = 1500cmV"^s"\ 
The values of carrier concentrations used are depicted in the 
plot. The spin-relaxation time, Ts, is set to 10~^s. The source 
is set at a distance of 2 x 10~^cm from the Fe/GaAs interface 
and it has a width of 5 x 10"'* cm. The source function is 
constant within this interval. 



FIG. 5: (Color online) The calculated spin detection sen- 
sitivity for an Fe/MgO/GaAs electrical spin detector with 
various levels of doping. On both panels, the red/gray solid 
line with circles is the Fe/MgO tunneling current spin polar- 
ization calculated from first-principles. The left panel is for 
/i = SOOOcmV'^s"^ and the right for /i = 1500cm^V~^ . 
The values of carrier concentrations used are depicted in the 
plot. The spin-relaxation time, Ts, is set to 10~^s. The source 
is set at a distance of 2 x 10~^cm from the Fe/GaAs interface 
and it has a width of 5 x 10~*cm. The source function is 
constant within this interval. 



semiconductor is more lightly doped. More specifically, 
from equations ((43|l and ([37]) we can see that the rate of 
drop/increase depends on the conducting properties of 
the semiconductor, the spin polarization of the interface 
(prefactor (1/G| — l/Gx)) and the location and width 
of the spin source (prefactor F{j)). The rate of increase 
during spin collection is proportional to cx n~^^^ while 
the rate of decrease is proportional to cx ji^^n^^^^ 

In Fig. [3] we show the influence of the tunneling resis- 
tance. In a ferromagnet /semiconductor interface the tun- 
neling resistance depends on the Schottky barrier height, 
width and shape. The first is determined by the magni- 
tude of the band gap in the semiconductor and the posi- 
tion of the Fermi level relative to the top of the valence 
band while the last two are modulated with doping. In 
Fig. [3] we can see that smaller tunneling resistance results 
in bigger variation of spin detection sensitivity with bias. 
This can be understood by from the linearized results in 
Eq. (1121). 

The rate of change of spin detection sensitivity with 
bias is inversely proportional to the tunneling resistance. 
The physical origin of this lies in that, for a given 
applied voltage, larger tunneling resistance will result 
in smaller current; this is reducing the eifect of drift. 
Generally speaking, the current spin polarization and 
hence (l/Gj — l/G^) can have a strong bias dependence. 
It was shown in Ref. '^.Ts'.l? that the current spin po- 
larization of Fe/GaAs(001) junctions has a very strong 



bias dependence, and even reverses sign within a small 
interval around zero bias. In Ref. [l9,20 two different 
microscopic models to explain the experimentally ob- 
served bias dependence of the tunneling current were dis- 
cussed. The bias dependence of spin sensitivity is due to a 
combination of the macroscopic physics described above 
and the microscopic bias dependence of (1/G| — l/Gj). 
To predict the resulting behavior in specific ferromag- 
net/semiconductor junctions we have incorporated our 
model first-principle results for the bias dependence of 
(1/GT-1/Gx)(y). 

In Fig. 2] we show the calculated spin detection sen- 
sitivity and first-principles current polarization for a 
Fe/GaAs(001) interface. The first-principle method and 
results are identical to that presented in Ref. [l^. As it 
was explained in Ref. [l^ the interface electronic struc- 
ture of Fe/GaAs(001) results in two strong minority-spin 
peaks in the energy dependence of electron transmission 
across the interface. One is located at about 125 meV 
below Ep and the other 125 meV above. The result is 
a strong bias dependence of the current spin polariza- 
tion. We see in Fig. [3] that the bias dependence of spin 
detection sensitivity bares some resemblance to that of 
the current spin polarization but in the general case can 
be significantly different. The magnitude of dV/dC de- 
creases (increases) faster than the spin polarization in 
the negative (positive) bias and the difference between 
the two increases as we make the semiconductor less con- 
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ductive. As we can see in Fig. |4]the spin detection sen- 
sitivity can be raised by an order of magnitude in the 
positive bias. Therefore, under certain conditions the 
macroscopic factors described above can have a domi- 
nant influence over the microscopic factors that influence 
the bias dependence of current polarization. Since it is 
much easier to control the conducting properties of semi- 
conductor rather the microscopic electronic properties of 
the interface, these effect can have a direct application in 
optimization of electrical spin detectors. This result is in 
very good agreement with the experimental bias depen- 
dence of spin-detection sensitivity presented in Ref. [la- 
It is interesting to examine the bias dependence of 
spin-detection sensitivity for a different interface. The 
Fe/MgO(001) interface is different from Fe/GaAs(001) 
in many ways and it is of great interest to spintronics 
community. To calculate the spin-dependent tunneling 
conductance of Fc/MgO(001) we used the same approach 
as in Ref. 19,21,22 and the same setup with Ref. [23. The 
details of the calculation such as the chosen k-mesh in the 
two dimensional Brillouin zone (2DBZ) and the method 
of calculation for the total current and spin polarization 
are the same with those used for Fe/GaAs(001) inter- 
face and are described in Refs. mH It was shown in 
Ref. [23112^ that the interface minority-spin resonances in 
Fe/MgO(001) interface are located far from the F point 
in the 2DBZ, contributing much less to the tunneling 
conductance than they do in the case of Fe/GaAs(001) 
interface. Therefore as we can see in Fig. [S]the spin po- 
larization of the tunneling current has less dramatic bias 
dependence in this case. The band gap of MgO is about 5 
times larger than the band gap of GaAs so the Fe/MgO 
tunneling barrier has much bigger resistance than the 
Fe/GaAs. Because of these differences with Fe/GaAs 
interface the spin-detection sensitivity is much less sensi- 
tive to the changes in applied bias. This is consistent with 
the analysis given so far and in particular with Eq. (|43p . 



IV. CONCLUSION 



of Fe/GaAs Schottky tunnel contacts show a very good 
agreement with experimenti^. As in the experiment we 
were able to enhance the spin sensitivity by an order of 
magnitude when applied positive voltage. Our results 
for Fe/MgO/GaAs show a similar enhancement though 
the magnitude of the effect is smaller than in Fe/GaAs. 
This is explained by the bigger height of tunneling barrier 
in the case of Fe/GaAs. These results suggest specifled 
routes on how to engineer efficient electrical spin detec- 
tors using ferromagnet/semiconductor interfaces. 
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V. APPENDIX: VOLTAGE DERIVATIVE OF 
THE INTERFACE CONDUCTANCES 

In the previous sections we neglected terms propor- 
tional to the voltage derivative of the interface conduc- 
tances because these terms are small for typical parame- 
ter values. In this appendix we discuss the contribution 
of these terms. For notational simplicity, it is convenient 
to define 

and 

In this notation the voltage drop at the interface is 



We presented a theory of electrical spin detection in 
ferromagnet/semiconductor detectors. We showed that 
the sensitivity of such detectors can have a strong bias 
dependence. The origin of this dependence lies in the 
microscopic electronic structure of the interface and the 
macroscopic electrical properties of the conducting chan- 
nel in the semiconductor. The first was incorporated 
in our model with the help of first principles electronic 
structure calculations. With the help of a model spin de- 
tector which has constant current polarization with re- 
spect to bias we showed that the latter by itself is ca- 
pable of producing strong bias dependence of sensitivity. 
This result suggests that enhancement of detector's sen- 
sitivity is possible independent of the materials used to 
construct the detector by engineering the electrical prop- 
erties of the conducting channel in the semiconductor 
and tuning the bias. Our results for the particular case 



R-P) 



(46) 



where P is the current density spin polarization P = 
(2/3o — 1). Including the voltage derivative of the interface 
conductances gives the voltage signal as 



dS^ 4 



diR++R-P) \ dy_ dP_ 
dV ) dS^ ^ dSb 



(47) 



where y ^^^^-^ ^^ j ^ is the new term. Solving for ^ 
gives 



dV_ 

dSb 



1^- dS'' 



j_ f d{R++R^P) \ 
4 \ dV J 



(48) 
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Equations ^ and ^ are used to find j^d(R+P+R^ 



3 dP V 



-, _ j d(R++R-P) ■ (^1) 
-^4 dV 



where 



2Ef 
~3e" 



A = i?. 



(1 + 250)1/3 (l-2(5")l/3 



2A_ 



^i^O') (50) 



(l + 2(5")i/3 (1-2(50)1/3 



Equations (gB]), (gH]), ([501), and dHI) give the voltage sig- 
nal including the voltage derivative of the interface con- 
ductances. The new terms that contain ''■(^■++^-^^ 

d(it'.+ P+-R-) 



dV 



dV 

sidered here 



are numerically small for the parameters con- 
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